Observation of the Shadowing of Cosmic Rays by the Moon using a Deep 

Underground Detector 



The MACRO Collaboration 

M. Ambrosiol2, R. Antolini^, C. Aramo^'", G. Auriemmai'*''', A. Baldinii^, G. C. Barbarinol^, B. C. Barish", G. Battistoni'^''', 
R. Bellotti\ C. Bemporad , P. Bernardini 1°, H. Bilokon ^, V. Bisi C. Bloise ^, C. Bower S. Bussino^"', F. Cafagna\ 
M. Calicchio\ D. Campana , M. Carbonic, M. Castellanoi, S. Cecchmi2><:, F. Cei i^'l^^ y. Chiarellaf*, B. C. Choudhary-*, S. Coutu 
L. De Benedictisl, G. De Cataldoi, H. Dekhissi ^.17^ q Marzo\ I. De Mitri^, J. Derkaoui ^'^'^ , M. De Vincenzi^^.e^ 
A. Di Credico , O. Erriquez^, C. Favuzzi^, C. Forti^, P. Fusco^^, G. Giacomelli ^, G. Gianninil^'-'', N. Giglietto^, M. Giorgini 2, 
M. Grassi L. Gray ^'^j A. Grillo , F. Guarino P. Guarnaccia^ C. Gustavino , A. Habig 3, K. Hanson R. Heinz Y. Huang'', 

E. larocci E. Katsavounidis* , E. Kearns ^, H. Kim^, S. Kyriazopoulou* , E. Lamanna C. Lane ^, D. S. Levin P. Lipari 
N. P. Longley M. J. Longo F. Maaroufi 2. 17^ G. Mancarella 1'', G. Mandrioli , S. Manzoor 2.™, A. Margiotta Neri 2, A. Marini 
- . , D. Martello A. Marzari-Chiesa M. N. Mazziotta^, C. Mazzotta 1°, D. G. Michael*, S. Mikheyev ■*>7,'i^ l. Miller P. Monacelli ^, 
OO I T. Montarulii, M. Monteno S. Mufson J. Musser D. Nicolol^.d, C. Orth 3, G. Osteria M. Ouchrif 2.17, Q. Palamara 

G^ V. Patera '^'S , L. Patrizii 2, R. Pazzi C. W. Peck", S. Petrera ^, P. Pistilli V. Popa V. Pugliese A. Raino^, J. Reynoldson 
Q\ 7, p Ronga 6, U. Rubizzo c. Satriano L. Satta ^-9, E. Scapparone 7, K. Scholberg 3, A. Sciubba "^'f , P. Serra-Lugaresi 2, 

M. Severi l*, M. Sioli 2, M. Sitta p. Spinellil, M. Spinetti "5, M. Spurio 2, R. Steinberg^, J. L. Stone 3, L. R. Sulak 3, A. Surdo 1°, 
G. Tarlel\ V. Togo 2, D. Ugolotti 2, M. Vakili C. W. Walter 3, and R. Webb l^, 

1. Dipartimento di Fisica dell'Universita di Bari and INFN, 70126 Bari, Italy 
^ ' 2. Dipartimento di Fisica dell'Universita di Bologna and INFN, 40126 Bologna, Italy 

, 3. Physics Department, Boston University, Boston, MA 02215, USA 

, 4. California Institute of Technology, Pasadena, CA 91125, USA 

5. Department of Physics, Drexel University, Philadelphia, PA 19104, USA 

6. Laboratori Nazionali di Frascati dell'INFN, 00044 Frascati (Roma), Italy 
7. Laboratori Nazionali del Gran Sasso dellTNFN, 67010 Assergi (L'Aquila), Italy 

8. Depts. of Physics and of Astronomy, Indiana University, Bloomington, IN 47405, USA 
9. Dipartimento di Fisica dell'Universita dell'Aquila and INFN, 67100 L'Aquila, Italy 
10. Dipartimento di Fisica dell'Universita di Lecce and INFN, 73100 Lecce, Italy 
^ . 11. Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 

^ ' 12. Dipartimento di Fisica dell'Universita di Napoli and INFN, 80125 Napoli, Italy 

I ' 13. Dipartimento di Fisica dell'Universita di Pisa and INFN, 56010 Pisa, Italy 

^ I 14. Dipartimento di Fisica dell'Universita di Roma "La Sapienza" and INFN, 00185 Roma, Italy 

D I 15. Physics Department, Texas A&M University, College Station, TX 77843, USA 

(~| , 16. Dipartimento di Fisica Sperimentale dell'Universita di Torino and INFN, 10125 Torino, Italy 

17. Also Faculty of Sciences, University Mohamed I, B.P. 424 Oujda, Morocco 
a Also Universita della Basilicata, 85100 Potenza, Italy 
' b Also INFN Milano, 20133 Milano, Italy 

rN ', c Also Istituto TESRE/CNR, 40129 Bologna, Italy 

d Also Scuola Normale Superiore di Pisa, 56010 Pisa, Italy 
e Also Dipartimento di Fisica, Universita di Roma Tre, Roma, Italy 
/ Also Universita di Trieste and INFN, 34100 Trieste, Italy 
g Also Dipartimento di Energetica, Universita di Roma, 00185 Roma, Italy 
h Also Institute for Nuclear Research, Russian Academy of Science, 117312 Moscow, Russia 
i Also Institute for Space Sciences, 76900 Bucharest, Romania 
/ Swarthmore College, Swarthmore, PA 19081, USA 
m RPD, PINSTECH, P.O. Nilore, Islamabad, Pakistan 
n Also INFN Catania, 95129 Catania, Italy 
o Also Department of Physics, Pennsylvania State University, University Park, PA 16801, USA 



o 
o 

OO 

0^ 



Using data collected by the MACRO experiment during the years 1989-1996, we show evidence 
for the shadow of the moon in the underground cosmic ray flux with a significance of 3.6cr. This 
detection of the shadowing effect is the first by an underground detector. A maximum-likelihood 
analysis is used to determine that the angular resolution of the apparatus is 0.9° ± 0.3°. These 
results demonstrate MACRO'S capabilities as a muon telescope by confirming its absolute pointing 
ability and quantifying its angular resolution. 
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I. INTRODUCTION 



MACRO is a large area underground detector located in Hall B of the Gran Sasso National Laboratory (LNGS) 
in Italy at an average depth of 3700 m.w.e. The full apparatus has dimensions 76.5mx 12nix9.6m. The detector's 
active technologies include liquid scintillation counters and streamer tubes. Low-radiation crushed rock separates the 
active detector planes. Normally incident muons from above lose « 1.5 GeV in traversing the detector; at the surface, 
throughgoing muons have energies in excess of ~ 3 TeV, corresponding to primary cosmic rays of energy > 7 TeV. 
The average energy of the muons detected by MACRO is ~ 300 GeV. A detailed description of the detector can be 
found in 

MACRO was primarily designed to search for monopoles and rare particles in the cosmic rays, including high energy 
neutrinos and muons from cosmic point sources 0] . Since high energy neutrinos are identified as upward-going muon 
secondaries from neutrino interactions in the rock below the detector, MACRO functions as a muon telescope in its 
search for muons and neutrinos from cosmic sources. The purpose of this investigation is to verify its absolute pointing 
and quantify its angular resolution to point sources. Confidence in pointing is clearly essential for a telescope, and a 
preliminary determination of MACRO'S pointing has been made by bootstrapping highway surveys from the surface 
into the tunnel. However, a check of the accuracy of this determination is important. Further, reconstructed muon 
tracks point back to an area on the sky whose width depends on the intrinsic angular resolution of the detector as well 
as multiple Coulomb scattering in the rock overburden. Consequently, the quality of MACRO as a muon telescope 
depends on its angular resolution. As the required search region around a source shrinks with improving resolution, 
the contribution of events from the flat or slowly varying background also decreases relative to the signal. 

Traditional astronomical telescopes use observable point sources to determine pointing and angular resolution. 
Despite early reports of their existence j^, however, there are as yet no established cosmic muon sources. In the 
absence of cosmic sources, we adopt the approach of the CYGNUS g, CAS A [|, and Tibet g air shower arrays who 
used the moon as a fiducial object. By blocking cosmic ray primaries, the moon appears as a cosmic ray antisource 
or "shadow" g which can be used to verify MACRO'S absolute pointing and to determine its angular resolution. 

II. THE MUON DATA SAMPLE AND THE EXPECTED BACKGROUND 

A. Muon data sample 

MACRO consists of six nearly identical units called supermodules, each of dimension 12.6mx 12mx9.6m. The lower 
4.8m of each supermodule consists of 10 horizontal planes of streamer tubes of dimension 12mx6m. The 8 innermost 
planes are separated by seven layers of 60 gm/cm^ absorber. The two outermost planes are each separated from the 
next nearest streamer tube plane by a 25 cm layer of liquid scintillator. The lateral walls consist of stacked tanks 
of liquid scintillator, 25 cm thick, sandwiched between six vertical streamer tubes planes. The upper 4.8 m of each 
supermodule (the attico) is a hollowed out version of the lower supermodule. It consists of a top plane of liquid 
scintillator with two streamer tube planes on both above and below, and lateral walls of liquid scintillator with three 
vertical streamer tubes planes outside and inside. All streamer tube wires are read out, providing the X coordinate 
on the horizontal planes and the Z coordinate on the vertical planes. On the horizontal planes the second coordinate, 
D, is obtained by horizontal aluminum strips oriented 26.5° with respect to the streamer tube axis. These strips allow 
stereoscopic reconstruction. For a complete description of the detector, see 

The muon sample used for the present analysis includes all events collected from the start of MACRO data taking in 
February, 1989 through the end of 1996. The sample totals 39.3 x 10^ events collected over 2.7 x 10"^ live days. During 
the first part of this period the apparatus was under construction. Long running periods included one supermodule 
without the attico {A^ff^l « 1,010 m^ sr), six supermodules without attico {Ae/ffl « 5,600 m^ sr) and finally the 
full six supermodules with attico {Aefffl « 6, 600 m^ sr). Approximately 60% of the data sample was obtained during 
periods when MACRO had full acceptance. 

The criteria used to select events for this analysis were designed to optimize the quality of reconstructed tracks. 
The selected events are consequently those which most accurately point back to their origin on the celestial sphere. 
The specific data cuts used in this analysis are listed below. 

• Run cuts 

1. Runs less than one hour in length were cut since short runs were usually abnormally terminated and often 
contained malfunctioning hardware or software. 

2. Runs with large numbers of UT clock errors were cut. 
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• Event cuts 



1. Only single and double muon events were retained; events with multiplicity > 2 were cut. The more 
complicated task of reconstructing multimuon tracks is more likely to introduce tracking errors. 

2. Successful track reconstruction requires a minimum of four crossed horizontal streamer tube planes; events 
with fewer than four crossed horizontal planes were cut. 

3. Events crossing fewer than 3 streamer tube planes in the lower supermodules were cut. This cut removes 
low energy events that pass through the attico without crossing any rock absorber layers. Such events have 
large multiple Coulomb scattering angles and so are spread widely with respect to their point of origin on 
the celestial sphere. 

4. Reconstructed tracks with y^/d.o.f. > 1.5 were cut. This cut removes events with poorly reconstruced 
tracks as well as events with large numbers of hits outside the track. 

5. Events with different reconstructed multiplicity in the two streamer tube views were cut. 

6. Events with UTC clock errors were cut since the correct time is necessary to project the track back onto 
the celestial sphere. 

These selection criteria reduce the sample size to 30.51 x 10^ muons. 

The topocentric position of the moon was computed at the arrival time of each event in the sample using the 
database of ephemerides available from the Jet Propulsion Laboratory, JPLEPH |Q. A correction for the parallax 
due to macro's instantaneous position on the earth was applied to each ephemeris position Since the parallax 
correction requires an accurate computation of the earth-moon distance, this calculation also results in an accurate 
determination of angular size subtended by the moon at MACRO. 

The muon events in a window 10° on a side and centered on the moon were selected for further analysis. There are 
2.3 X 10^ events that pass all cuts in this window. 



Twenty five background samples were generated for each run used in the analysis. These backgrounds were 
constructed by coupling the direction of each muon in the run with the times of 25 randomly selected muons from 
the same run. The 25 background samples were then processed using the same procedure as the muon data sample. 



Each muon event in the window was sorted into bins of equal angular width as a function of angular distance from the 
moon center. The angular width of each bin was 0.125°, which gives a solid angle for the ith bin of Afij = (2i— 1) -0.05 
deg^. Once filled, the contents of each bin, Ni, were divided by the solid area of the bin, resulting in the distribution, 
AiV/Af2 which approximates the differential event density as a function of angular distance from the moon center. 
Figure |l] shows as data points the computed differential event density. The errors on the data points are statistical. 
Superposed on this distribution is the expected event density in the absence of the shadowing effect, AN'''^'^ / AD,. 
This background distribution was determined by averaging over the 25 background samples the number of events 
in each solid angle bin, and then dividing the result by the solid angle of the bin. This nearly flat distribution is 
described reasonably well by a constant event density « 735 deg~^. The shadowing effect, or the deviation of the 
angular event density from the background, is clearly evident. This figure confirms that reconstructed muon tracks 
can be accurately pointed back to celestial coordinates, thus confirming MACRO'S absolute pointing ability. 

We have made a simple estimate of the significance of this detection of the moon shadow using the information 
given in Figure |^. First, we computed the integral event deficit as a function of the angular distance from the moon 
center. Out to bin n, the integral event deficit, AnN'^'^^ , is given by 



B. Expected background 



III. SHADOW OF THE MOON 



A. Event deficit around the moon 



n 




(3.1) 
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This distribution represents the cumulative number of events that the data distribution differs from the flat background 
and is shown in Figure ^ This figure shows that the deficit increases until the shadow can no longer be distinguished 
from the background at ~ 0.9°. Out to 0.9° the integral deficit can be approximated by /^N'^'^^ k, 165a, where a 
is the angular distance from the moon center; this simple approximation is shown as a dashed line in Figure ^. An 
estimate of the significance is then given by 

ATV'*'^/ 165a ^ ^ 

, = , « 3.5. (3.2) 

^/N^ V7357ra2 

This computation shows that the moon shadow detection has a significance of approximately 3.5(t. 



B. Maximum likelihood analysis 

In the simple deficit analysis above, we have implicitly assumed that the position of the moon's shadow is known 
and we have binned the events using this information. We now relax that assumption and search for the moon shadow 
in a direction-independent way with the maximum likelihood method of COS-B |l0[| , a technique first described in 
detail by Cash ||ll[] . This method is based on a priori knowledge of the point spread function of the MACRO detector 
(MPSF). 

We have determined the MPSF using the observed space angle distribution of double muons. As shown in [|l2|, the 
number of double muons as a function of lateral separation is a power law distribution that falls with a scale length 
of 15 m. Muon pairs produced in a primary interaction at 20 km therefore have typical initial separation angles 
< .05°. Since "double" muons initially move along virtually parallel paths, the distribution of their separation angles 
is a good measure of the deviations introduced into their tracking parameters by both scattering in the mountain 
overburden and the detector's intrinsic angular resolution. The space angle distribution of "double" muons must be 
divided by \/2 to obtain the MPSF since both muons deviate from their initial trajectories. In Figure ^ the MPSF, as 
determined from 1,044,877 muon pairs, is shown in altitude and azimuthal coordinates. Figured shows the strongly 
peaked, non-Gaussian behavior of the MPSF. 

To find the most likely position of the moon, we compare the two dimensional distribution of muons in the window 
centered on the moon with the expected background events in the same window. In this analysis, each muon event is 
first sorted into a grid of equal solid angle bins (Afi = 0.125° x 0.125° = 1.6 x 10~^ deg^). The shadowing source of 
strength Sm at fixed position (x^, ys) that best fits the data is then found by minimizing 

x\xs,ys, Sm) = 2 ^ [iV^ -N, + N, In -^], (3.3) 

where the sum is over all bins in the window [ p^ . Here Ni is the number of events observed in each bin j, N!^^ is 
the expected number of events in bin i, and nt,i„ is the number of bins in the grid. This expression assumes that a 
Poissonian process is responsible for the events seen in each bin. The expected number of events in bin i is given by 

Nr = N^"" - Sm ■ V{x, - x,,ys - yi), (3.4) 

where N^'^'^ is the average number of background events at position (xi, yi), and Sm • "Pixs — Xi,ys — yi) is the number 
of events removed from bin i by the shadow of the moon. Here V{xs — Xi,ys — yi) is the MPSF, modified for the 
finite size of the moon's disk, computed at the point {xi, yi) when the shadowing source is at {xs,ys). The MPSF was 
modified by first selecting random positions on a disk with the average lunar radius, 0.26°, and then drawing offsets 
from these positions from the MPSF distribution These new positions were rebinned into a new histogram and the 
resulting distribution normalized to unit area. The modified point spread function is shown in Fig. |^. In this figure, 
the unmodified, normalized MPSF is shown in Fig. the MPSF modified for the moon's finite disk is shown in 
Fig. In Figs. |c and d the MPSF has been modified by still larger shadowing disks. The effect of the finite size 
of the moon's disk does not have a large impact on the analysis. 

Finally, the shadow strength Sm that minimizes "^as computed for every grid point in the window. This minimum 
x'^{xs,ys, Sm) was then compared with x^(0) for the null hypothesis that no shadowing source is present in the window 
{Sm = 0), 

\ = x\0)-xHxs,ys,SM). (3.5) 

The most likely position of the moon is the bin in which the maximum A, A™"^ = A, is found. Since there is only 
one free parameter, Sm, A behaves like Xi, ^X^ distribution with one degree of freedom pO[ |. The significance of the 
moon detection is given by Xi(A)- 



4 



In Figure || we show the results of this analysis in a window 4.375° x 4.375° centered on the moon. This window 
has been divided into 35 x 35 cells, each having dimensions 0.125° x 0.125°. In this figure, A is displayed in grey 
scale format for every bin in the moon window. Also shown is the fiducial position of the moon and a circle centered 
at this position corresponding to the average lunar radius, 0.26°. The maximum A — 18.3 is found somewhat offset 
from the fiducial moon position at A Azimuth = —0.25°, A Altitude = +0.125° and provides further confirmation of 
macro's absolute pointing. The value of the shadow strength at this position, Sm = 153.9 ± 37 events, agrees well 
with the observed value 155 events. Although the displacement from the fiducial position of the moon is not of high 
significance, we note that it is consistent with the deflection of cosmic ray primaries by the geomagnetic field 

We have verified the properties of the A distribution by constructing 71 other windows similar to the moon window, 
each displaced from the next by 5° in right ascension. For each off-source window, we followed the procedure used 
for the moon window in computing the expected background. To avoid edge effects associated with a source near 
the edge of a window, we only evaluated A for the central 12x12 bins. In Figure ^ we have plotted the cumulative 
distribution of the 12 x 12 x 71 = 10, 224 values of A for the 71 off-source windows. Superposed on this distribution 
as a solid curve is the cumulative Xi distribution. The distributions agree well (the structure seen at 7 < A < 10 
is due to roundoff error). Thus, the probability of the detection of the shadow of the moon at this bin location is 
Pixi — 18.3) < 1.9 X 10^^. However, we would have considered the detection of the moon shadow equally secure had 
the maximum been found at any of the 9 bins within the geometric shadow, which reduces the probability of detection 
to p < 1.7 X 10~^. The significance of the detection is therefore 3.6a and the results of the likelihood analysis and the 
deficit analysis are in excellent agreement. 

As a crosscheck, a similar analyis in right ascension and declination was performed. The results were equivalent, 
as required. 



IV. ANGULAR RESOLUTION OF THE MACRO APPARATUS 

As shown in Figure ||, the MPSF cannot fit by a simple Gaussian function. Thus, the technique used by the 
CYGNUS 1^, CASA and Tibet ^ air shower arrays, where such a simplification was possible, must be modified. 
The air shower experiments find their PSF by rigidly scaling the dispersion of their Gaussian resolution functions 
and then computing the likelihood function for the detection of the moon shadow for each scaled value of a. The 
maximum likelihood defines the a to be used in the computation of the angular resolution. In our approach, we first 
defined a scale parameter T that rigidly scales the modified MPSF by the factor T , 

V{xs - Xi,ys - yi;T). 

We then repeated the likelihood analysis in the moon window for different values of T . We assume that the value of 
T that maximizes A gives the best V for computing the angular resolution. 

The distribution of A as a function of T is shown in Figure 0. The maximum A, A* is found ior T = 1.0 which 
implies that the unsealed space angle distribution of double muons should be used to determine the anglular resolution 
of the MACRO apparatus. Many definitions have been used for the angular resolution ( [0, ||], |^). We choose the 
cone of angle 6'gg% that contains the 68% of the events from a point-like source. In Figure]| we show the space angle 
distribution of double muons. Using our definition of the angular resolution. Figure gives 6*58% — 0.90°. The Icr 



error limits on the angular resolution can be estimated from the interval of T in which A* falls by 1.0 1 13|. The shaded 
region in Figure |7| shows this interval to be = 1.0 ± 0.35. Using these values to scale the double muon distribution, 
we find that the angular resolution is 6'g8% = 0.9° ± 0.3°. 

In Figure ^ we show as a solid line the expected event density in the moon window for an angular resolution of the 
MACRO apparatus of 0.90°. The model used in this computation is given by 

^ = A[Ar'"='^ - Sm ■ Via; T = 1.0)]/Af], (4.1) 



where a = \J (xs — Xi)^ + {ys — yi)^ is the angular distance from the moon center. The model fits the data well 
(X^ = 48.7/47 d.o.f.). 



V. CONCLUSIONS 



The MACRO detector at a mean depth of 3700 m.w.e, operational since February 1989, has collected a muon sample 
of about 39 million events. Using this sample we have searched for the moon shadow cast in the cosmic ray sky at 
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primary energies ^ 10 — 15 TcV. In the deficit analysis, we find an event deficit around the moon of significance 3.5(t. 
With a maximum likelihood analysis, we confirm the detection of the moon's shadow with significance of 3.6cr. This 
is the first detection of the moon shadow underground. Our estimate of the angular resolution is 6*68% = 0.9° ± 0.3°. 

These results demonstrate MACRO'S capabilities as a muon telescope by confirming its absolute pointing ability 
and quantifying its angular resolution. This investigation shows that the MACRO detector has the capability of 
detecting signals from cosmic sources by observing secondary cosmic muons underground. 
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IS 2 2.5 3 3.5 4 

Angular distance from the Moon (degrees) 

FIG. 1. The event density vs angular distance from the moon center in bins of equal angular width. The width of each bin is 
0.125°. The dashed curve is the average expected background computed from 25 background samples. The solid curve shows 
the expected event density as computed for an angular resolution of the MACRO apparatus of 0.90°. 
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FIG. 2. Integral event deficit vs angular distance from the moon center. Superposed on the integral event deficit is an 
approximate fit to the distribution out to 0.9°: AN"^"-^ = 165a, where a is the angular distance from the moon center. 



8 




FIG. 3. (a) The point spread function of the MACRO apparatus (MPSF) derived from the space angle distribution of 
1,044,877 double muons. (b) Contour diagram of the MPSF. The contours levels are equally spaced between 2,000 and 18,000 
muons. The MPSF is circularly symmetric in this coordinate system. 
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<c) (d) 

FIG. 4. (a) The MPSF, normalized to unit area from Fig. ||a. (b) The MPSF modified by a disk of radius 0.26°, the average 
lunar radius, and then normalized to unit area; this is the modified MPSF used in the likelihood analysis, (c) The MPSF 
modified by a disk of radius 0.4° and then normalized, (d) The MPSF modified by a disk of radius 1.0° and then normalized. 
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FIG. 5. The two dimensional distribution of A in bins of equal solid angle in the moon window The axes are offsets from 

the moon center. The fiducial position of the moon, at position (0,0), is marked by a +; a circle corresponding to the average 
lunar radius, 0.26°, is centered at this position. The A grey scale is given at the right margin of the figure. The maximum of 
this distribution, A = 18.3, is offset from the fiducial moon position at AAzimuth = —0.25° and A Altitude = +0.125°. 



11 




2 4 6 8 10 12 14 16 18 20 

^max 

FIG. 6. The cumulative distribution of 10,224 values of A for all bins in the 71 off-source windows. Superposed on this 
distribution, as a solid curve, is the cumulative Xi distribution. The arrow marks A = 18.3 from the likelihood analysis. 
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degrees 

FIG. 8. The space angle distribution of double muons. The angular resolution of the MACRO apparatus, defined as the cone 
angle containing 68% of the events from a point source, is shown as a shaded region. With this definition, MACRO'S angular 
resolution is 0.90°. 
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